Spin-dependent forces on trapped ions for phase-stable quantum gates and motional 

Schrodinger-cat states 
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Favored schemes for trapped-ion quantum logic gates use bichromatic laser fields to couple internal 
qubit states with external motion through a "spin-dependent force." We introduce a new degree 
of freedom in this coupling that reduces its sensitivity to phase decoherence. We demonstrate 
bichromatic spin-dependent forces on a single trapped 111 Cd + ion, and show that phase coherence 
of the resulting "Schrodinger-cat" states of motion depends critically upon the spectral arrangement 
of the optical fields. This applies directly to the operation of entangling gates on multiple ions. 

PACS numbers: 42.50.Vk,03.67.Mn,03.67.Lx 



Trapped atomic ions have a number of desirable fea- 
tures that make them well suited for quantum informa- 
tion applications Q. Pairs of hyperfine ground states 
provide an ideal host for quantum bits (qubits) that can 
be manipulated and entangled via optical Raman transi- 
tions J2| ■ Since the first proposal for entangling trapped 
ions |3| , significant theoretical advances have relaxed the 
constraints on implementation to realize more robust en- 
tanglement schemes 0, 0, • These schemes rely on 
a spin-dependent force acting on each ion where "spin" 
refers to the effective Pauli spin associated with the ion- 
qubit's two- level system. Acting on a single ion, the force 
can entangle the ion's spin degree of freedom with its mo- 
tion; following disentanglement, the ion acquires a net 
geometric phase that is spin dependent @. When the 
force is applied to two ions, the geometric phase depends 
nonlincarly on their spins through their mutual Coulomb 
interaction, generating entanglement. 

Optical Raman fields are a convenient way to cre- 
ate strong spin-dependent forces for hyperfine ion-qubits 
For example, an entangling gate based on a b z - 
dependent force has been realized with a state-dependent 
AC Stark shift from Raman fields 0]j where b x ,y,z are 
the Pauli operators. Unfortunately, such b z gates are not 
compatible with magnetic field insensitive (or "clock") 
qubit states 0, 0, and are therefore open to qubit 
phase decoherence from fluctuating magnetic fields. An 
alternative solution is the M0lmer-S0rensen (MS) gate 
0, 0, , which uses a more complicated arrangement 
of bichromatic fields to realize a cr^-type force where b^ 
is a linear combination of b x and b y operators. While the 
gate does work on clock states, it can have a significant 
phase instability due to the hyperfine coherences required 
to generate the spin dependence. All gates, whatever the 
spin dependence, are susceptible to fast phase fluctua- 
tions of the Raman fields during the course of gate evo- 
lution. However, the MS gate can also be susceptible to 
slow phase drifts, which cause changes in the spin de- 
pendence of the force and results in a phase instability 
between gates. 

In this Letter we show how to overcome phase drifts in 



the MS gate, allowing clock-state benefits to be realized 
for ion-based qubits. Such ideas may be useful in other 
experimental settings since the MS scheme is a general 
paradigm for entanglement [jjj- A bichromatic force is 
first implemented on a single ion to demonstrate the pe- 
riodic formation of "Schrodinger cat" states of entangled 
spin and motion ||. The generation of single-ion cat 
states is a useful test-bed for issues related to the MS 
gate without the added complication of two ions. In par- 
ticular, we make use of the cat signal to demonstrate a 
phase-insensitive gate arrangement. 

A b(p force operates in a spin basis dressed by the Ra- 
man laser fields. It is created from a bichromatic field 
composed of Raman couplings to the first vibrational 
sidebands of the trapped ion's motion, assumed along the 
z-axis. The red sideband coupling | |,n + l) | T,n), 
and the blue coupling | I, n) <-> | | , create a Jaynes- 

Cummings and anti-Jaynes-Cummings interaction be- 
tween the ion's spin and its vibrational levels 
If the red and blue couplings have the same base Rabi 
frequencies (£l r = Of, = £l s b for n = 0) and balanced de- 
tunings (— S r = 5b = 6) then a spin-dependent, oscillating 
force is realized under the rotatin g-w ave and Lamb-Dicke 
approximations 0] (see also ref. The Hamiltonian 

in the interaction frame is 
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where z a — y h/2mui z is the size of the harmonic os- 
cillator ground state for an ion of mass m and center- 
of-mass oscillation frequency lo z . The strength of the 
force is given by F z = ftTL s b and and a are the os- 
cillator raising and lowering operators. The orientation 
of the force's spin dependence, b^ = er l ^ s b + + e^ s <7-_, 
is defined by the azimuthal angle <p s where b± are the 
spin raising and lowering operators. Note b^ = b x for 
the case 4> s = 0. In terms of the phases <p r and <pb of 
the red and blue sideband fields, the spin phase of the 
MS force in Eqn. ^ is 4> s = (0& + ^V)/2 and the mo- 
tional phase is <j) m — (4>b — 4'r)l^-- The unitary evolution 
generated by the force is a spin-dependent displacement 
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U{t) = V(a{tm^){U\+V{-a{t))\U)(U\ where 
and H0) are the eigenstates of 01 • IS the dis- 

placement operator in position-momentum (z-p) phase 
space with a(t) — cv e~^ m (l — e~ lSt ) and a — fl s b/2S. 
The corresponding position and momentum displace- 
ments are 2z Q Re(a(t)) and 2mu z z Im(a(t)) respectively. 

Instabilities in the spin and motional phases <f) s and <f) m 
arise from optical phase fluctuations 8(f> between the laser 
beams generating the Raman field at the ion [Fig. [I], 
or equivalently from fluctuations in the ion center-of- 
mass position. Slow global fluctuations in the force's 
motional phase are unimportant in the context of MS 
entangling gates, which involve only transient motional 
coherences associated with closed trajectories in phase 
space. Motional coherence between gates is therefore ir- 
relevant. Such is not the case for spin coherences; how- 
ever, a bichromatic field offers an additional degree of 
freedom to suppress the effect of fluctuations on (f> s . The 
two fields can be composed of co- or counter- propagat- 
ing Raman beat waves. For the co-propagating (phase 
sensitive) setup in Fig.^a), the phase fluctuation 5(f) ap- 
pears on the red and blue sideband couplings with the 
same sign ((f) r — (f>b = 5(f)) and so is directly written onto 
the MS spin phase ((f> s = 5(f)) . In the counter-propagating 
(phase insensitive) case in Fig.^b), the phase fluctuation 
appears with the opposite sign (— (f> r = (f>b — 6(f)) and is 
thus eliminated from the spin phase (</> s = 0). The phase- 
sensitive setup can be useful for cancelling common-mode 
fluctuations with other gate operations employing the 
same Raman beams while the phase-insensitive setup is 
useful for synchronizing with other operations involving 
microwave fields or different Raman beam geometries. 
This discussion only accounts for fluctuations shared by 
the red and blue sideband fields as shown in Fig. ^ but 
independent fluctuations can also be avoided with an ap- 
propriate optical design. 

Our experiment is centered around a linear ion-trap 
[lij employing both rf and dc fields to create a harmonic 
confining potential, characterized by trapping frequen- 
cies {u) x ,oj y ,uj z }/2ir = {8,9,3.55}MHz. Single ni Cd+ 
ions are loaded into the trap from a neutral CdO oven 
source via resonant pulsed-laser photoionization. The 
qubit resides in the 5Si/ 2 ground hyperfine clock states, 
||) = \F = l,m F = 0) and |f) = \F = 0,m F =0), sep- 
arated by u)i l f/2n = I4.53GHz. Doppler cooling on the 
5Si/2 — 5P3/2 cycling transition (wavelength 214. 5nm, 
linewidth j/2ir — 47MHz) brings the ion to a thermal 
state given by n w 6. An optional step of sideband cool- 
ing can further cool the ion to near the ground state 
of motion along the £-axis (n « 0.05) 0, well within 
the Lamb-Dicke limit. Optical pumping allows near per- 
fect preparation of the qubit to the state |T). High ef- 
ficiency qubit detection is accomplished through fluores- 
cence detection, where strong ion fluorescence is indica- 
tive of state |t) whereas state |J.) remains nearly dark. 

Raman transitions are generated with a pair of laser 
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FIG. I: Two possible Raman setups (top) to realize a bichro- 
matic force on a hyperfine ion-qubit and corresponding optical 
spectrum (bottom) showing red (rsb) and blue (bsb) sideband 
beatnotes. Two laser beams (solid lines, top) with optical car- 
rier frequency v generate (a) co-propagating or (b) counter- 
propagating Raman running waves (dotted lines) at the ion 
with wavevectors Ak ri b. Symbols include optical phase fluc- 
tuation 8(f>, hyperfine splitting u^f, ion trap frequency lu z , 
and force detuning 8. 



beams in a 90° geometry with wavevector difference Ak 
along the z-axis [Fig. Coupling occurs via the P 3 / 2 
state from which the beams are detuned by ~220GHz. 
A combination of laser modulation techniques is used 
to generate the Raman beat notes 0|. To bridge the 
hyperfine gap an electro-optic modulator is used to im- 
part a microwave frequency comb onto the Raman laser, 
which is then split into two arms. Since the comb exists 
on both arms, all frequencies exist to access beat waves 
propagating in either direction. A pair of acousto-optic 
modulators (AOMs), one in each arm, generates a fre- 
quency offset between the two beams. This allows for 
tuning across the various motional sideband resonances 
associated with either wavevector direction. A bichro- 
matic Raman field is obtained by driving one AOM with 
two rf frequencies [Fig.Q] so that any optical phase fluctu- 
ations on the red and blue sidebands are common mode. 
The sideband Rabi frequencies, typically 2kHz, are bal- 
anced in situ to better than 10% and the detunings 5 r 
and 5b to better than 100Hz. 

The action of the bichromatic force is demonstrated by 
applying it to a single ion, initially prepared in ||). This 
initial state can be expressed in terms of the diagonal 
spin basis of the bichromatic force as -^(If^.) + ||a;))|0), 
where without loss of generality the force is taken to be 
<T x -dependent and the initial motional state is the ground 
state. Application of the force generates the entangled 
state (|Tx) l a ) + |J.x) I - a ))i which can be written in 
the measurement or Jr, basis as |f) (\a) + |— a)) /2 + 
ID (l a ) — \~ a )) /2 Him Taking into account an initial 
thermal motional state as well as weak decoherence from 
motional heating [2l| leads to the general expression for 
the probability of measuring the ion to be in state ||): 



P, c (r) = -[l-e- 



tfi,r|4a | 2 -(n+i)|2a(T)| 2 
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FIG. 2: Single-ion evolution due to a spin-dependent Dichro- 
matic force, (a) Probability P? of measuring state \[) plot- 
ted vs. force duration r. Ion is initially Doppler cooled. 
Data is run-time averaged (100 shots/point) and smoothed. 
Solid line is a fit to Eqn. D2J modified to include a lin- 
ear change in peak and contrast (from spontaneous emis- 
sion) and a detuning drift across the data. Fit parame- 
ters include f2 si) /27r=2.2kHz(nxed), <5/27r=-5.46(3)kHz and 
n=8.1(3). (b)&(c)Phase space sketches of the ion motion at 
points indicated in (a). 



The initial thermal distribution is characterized by aver- 
age occupation number n, the motional heating rate by 
n and the duration of the force by r. 

An experiment in which the duration of the force is 
varied before qubit detection demonstrates the periodic 
entanglement in single ion evolution [Fig. [2]. A transition 
to = 1/2 indicates when a Schrodinger cat is formed. 
At this point, the motional wavepackets of the cat state 
are sufficiently far separated that the spin interference is 
inhibited, yielding equal probability of |t) and ||) in the 
measurement basis. At fringe minima corresponding to 
5t = 2rmr with m an integer, the motional states return 
to their original position and overlap. The spin inter- 
ference is restored, giving the initial state ||) such that 
?! k 0. This coherent process of periodic entanglement 
and disentanglement of the spin and motional degrees of 
freedom continues with reasonable signal quality for at 
least two oscillations [Fig.0- 

A complementary experiment in which the duration of 
the force is fixed but the detuning 6 is varied is shown 
in Fig [3J We have performed this experiment for two 
different initial temperatures of the ion, a "hot" case of 
a Doppler cooled ion and a "cold" case of a ion cooled 
to near the ground state. On resonance where the force 
is strongest, the inferred cat state separation in Fig. |3| 
is Az ss 10z o , a factor of 2.8 larger than the rms size 
of the "hot" ion's thermal state. The hot case has a 
broader envelope and narrower fringes. This occurs be- 
cause the average over the initial thermal distribution 
quickly draws the experiment outcome towards P£ = 1/2, 
even for small displacements. Nevertheless, all initial 
states should return on themselves at the same moment 
(within the Lamb-Dicke approximation) giving a full re- 
vival. The overall decrease of contrast, particularly vis- 
ible in the non-zero baseline, is due to spin decoher- 
ence and optical pumping induced by spontaneous emis- 
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FIG. 3: Probability P? plotted vs. detuning 5 of the bichro- 
matic force for (a) a ground-state cooled and (b) a Doppler 
cooled ion, initially prepared in The force is applied 

for 500^is. Data is run-time averaged with 100 shots/point 
and smoothed. Solid lines show fits to Eqn. |5] modified to 
include overall peak and contrast factors (for spontaneous 
emission) and a detuning drift across the data. An initial fit 
for data in (a) with n=0.05 fixed gives f2 IJ f,/27r=1.62(3)kHz 
and n=0.44(2)ms _1 . A subsequent fit for data in (b) with 
Q, s b/2-K— 1.62kHz fixed gives n=5.6(l) and n=0.62(6)ms~ . 
The values of n are 2-3 times larger than the directly mea- 
sured heating rate 0.2ms _1 |lq|. Phase space sketches (c)-(e) 
indicate ion evolution at detunings referenced in (a). 



sion. The detuning-dependent fringe contrast is consis- 
tent with motional decohcrcncc which has a character- 
istic exponential sensitivity to the size of the cat state, 
largest near resonance (6 ~ 0) [21I Eif . 

In order to observe the phase sensitivity of the MS 
force an interference technique is required. Ramsey Ln- 
tcrferomctry combined with Schrodinger cat formation 
provides a fluorescence signal that is sensitive to the ori- 
entation of the force's spin basis, characterized by the 
phase 4> s . The essence of the procedure is to rotate the 
spin of the initial state ||) with a 7r/2 pulse into the "x- 
y" plane at some azimuthal angle <p a . This results in the 
state l/\/2(|T) + e l ^°\i)) where <p is the reference phase 
associated with the 7r/2 pulse. The MS force is then ap- 
plied. Finally a 37r/2 analysis pulse returns the spin to 
the "z-axis" before measurement. A signal sensitive to 
d> s is obtained as follows: 



A=Pf(r) sin 2 



.) 



(3) 



where Pf ( r ) is given by Eqn. As long as the detun- 
ing and duration t of the force are chosen to generate a 
significant displacement (a ^> 1), the signal is approxi- 
mately ^sin 2 ((f> — <j> s ). If the initial 7r/2 pulse rotates 
the ion's spin into a state in which the force is diagonal 
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FIG. 4: (a) Interferometric photon-echo sequence to test op- 
tical phase sensitivity of the MS force. Phase <f> and duration 
of spin-rotation pulses (unshaded) indicated. The MS pulse 
(with spin phase <f> a ) and pulse for ac Stark shift compen- 
sation are shown shaded. (b)&(c) Probability P± plotted vs. 
applied shift in (j> where (b) and (c) correspond to setups in 
Figs. H^a)&(b) respectively. Data is run-time averaged with 
100 shots/point requiring ~200ms/point. Parameters include 
7r/2 pulse time of 13/iS, O s f,/27r«2kHz, 5/27r~5kHz, r=90^is, 
and n«6. Solid lines are a sinusoidal fit (b) and linear fit (c). 

{(f) = 4> s +TO7T with m an integer), a displacement occurs 
but no entangled cat state is formed. Therefore follow- 
ing the analysis pulse, the ion's spin returns to its initial 
state |t) such that P| = 0. On the other hand, if the 
rotated initial state deviates from this special condition 
{4>o 7^ 4>s + rmr), then the state is a superposition of the 
force's spin basis and a Schrodinger cat is formed. This 
is revealed by a net fluorescence signal on analysis. 

For convenience, we make use of an optical Raman car- 
rier transition to drive the tt/2 pulses. The actual pulse 
sequence for the interferometry employs a photon echo 
scheme where the the 3tt/2 analysis pulse is divided into 
a 7r and tt/2 pulse [Fig.^Ja)]. This provides a convenient 
way to cancel the effect of Stark shifts (~20kHz) from the 
optical carrier pulses as well as from the MS force itself. 
The Stark-induced phase from the MS pulse is cancelled 
by introducing a bichromatic pulse, far detuned from the 
motional resonance, into the second echo zone [Fig.^Ja)]. 

The carrier pulses which act as a phase reference are 
subject to the same optical phase drifts as the sidebands 
creating the bichromatic force. Therefore the interfero- 
metric signal Pj is expected to be stable for the phase 
sensitive Raman setup [Fig. ^a)] where fluctuations in 
the force's spin phase 4> s are common to the reference 
phase 4> . A scan of the reference phase [Fig.^Jb)] shows 
interference fringes with nearly full contrast (maximum 
being 1/2). This remains the case even with a piezoelec- 
tric transducer slowly (1Hz) modulating one Raman arm 
by an optical wavelength. Wash-out of the interferomet- 



ric signal should occur when the MS force is generated 
with the phase insensitive setup [Fig. ID[b)]. ^ n * ms case, 
the force's spin dependence, now decoupled from instabil- 
ities, is no longer correlated with the optical phase fluc- 
tuations on the carrier pulses. This configuration did not 
require the addition of any Raman path-length modula- 
tion to spoil the fringe contrast as inherent phase drifts 
on the optical table were a sufficient source of noise over 
the 200ms experiment averaging time. 

A spin-dependent force for clock-state qubits has been 
demonstrated that can be combined with other opera- 
tions in a way that is insensitive to slow optical phase 
fluctuations in the Raman beams and to drifts in ion 
center-of-mass position. This opens the way to phase- 
stable algorithms as follows: Single qubit rotations can be 
implemented with a co-propagating Raman transition or 
microwaves. The entangling gate may be realized directly 
with a er^-dependent force based on the phase-insensitive 
Raman setup. Alternatively, the phase-sensitive setup 
can be used in an echo sequence as in Fig. 01a) to con- 
struct a <7 Z gate piecewise where the carrier pulses act as 
a tr z <-> a x basis transformation. This arrangement can 
also suppress phase noise due to slow fluctuations in the 
relative position of two ions. 

This work was supported by the U.S. National Security 
Agency and Advanced Research and Development Activ- 
ity under Army Research Office contract DAAD19-01-1- 
0667 and the National Science Foundation Information 
Technology Research Program. 



[1] D. J. Wineland et al., J. Res. Nat. Inst. Stand. Tech. 103, 
259 (1998). 

[2] B. B. Blinov et al., Quant. Inf. Proc. 3, 1 (2004). 

[3] J. I. Cirac and P. Zoller, Phys. Rev. Lett. 74, 4091 (1995). 

[4] K. M0lmer and A. S0rensen, Phys. Rev. Lett. 82, 1835 

(1999) . 

[5] E. Solano, R. L. de Matos Filho, and N. Zagury, Phys. 

Rev. A 59, R2539 (1999). 
[6] G. J. Milburn, S. Schneider, and D. F. V. James, 

Fortschr. Physik (2000). 
[7] J. J. Garcia-Ripoll, P. Zoller, and J. I. Cirac, Phys. Rev. 

Lett. 91, 157901 (2003). 
[8] C. Monroe et al., Science 272, 1131 (1996). 
[9] D. Leibfried et al., Rev. Mod. Phys. 75, 281 (2003). 
[10] D. Leibfried et al., Nature 422, 412 (2003). 
[11] C. Langer and D. J. Wineland, private communication. 
[12] C. A. Sackett et al., Nature 404, 256 (2000). 
[13] A. S0rensen and K. M0lmer, Phys. Rev. A 62, 022311 

(2000) . 

[14] X. Wang, A. Sorensen, and K. Molmer, Phys. Rev. Lett. 

86, 3907 (2001); K. Helmerson and L. You, Phys. Rev. 

Lett. 87, 170402 (2001); S. Osnaghi et al, Phys. Rev. 

Lett. 87, 037902 (2001). 
[15] S. Wallentowitz and W. Vogel, Phys. Rev. Lett. 75, 2932 

(1995); S. Wallentowitz and W. Vogel, Phys. Rev. A 54, 

3322 (1996); S.-B. Zheng, Phys. Rev. A 58, 761 (1998); 



■5 



H. Moya-Cessa, S. Wallentowitz and W. Vogel, Phys. 
Rev. A 59, 2920 (1999). 

[16] E. Solano, G. S. Agarwal, and H. Walther, Phys. Rev. 
Lett. 90, 027903 (2003). 

[17] A global displacement-induced phase 9(t) has been omit- 
ted from U(t). 

[18] L. Deslauriers et al., Phys. Rev. A 70, 043408 (2004). 



[19] P. J. Lee et al., Opt. Lett. 28, 1582 (2003). 
[20] R. L. de Matos Filho and W. Vogel, Phys. Rev. Lett. 76, 
608 (1996). 

[21] Q. A. Turchette et al., Phys. Rev. A 61, 063418 (2000). 
[22] D. F. Walls and G. J. Milburn, Phys. Rev. A 31, 2403 
(1985). 



